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ABSTRACT 

We discuss the circumstances under which the bending of a jet can generate an internal 
shock wave. The analysis is carried out for relativistic and non-relativistic astrophys- 
ical jets. The calculations are done by the method of characteristics for the ca se of 



steady simple waves. This generalises the non-relativistic treatment first used by Ickc 



(1991). We show that it is possible to obtain an upper limit to the bending angle of a 
jet m order not to create a shock wave at the end of the curvature. This limiting angle 
has a value of ^ 75 for non-relativistic jets with a polytropic index k = 4/3, ~ 135 
for non-relativistic jets with k = 5/3 and ~ 50 for relativistic jets with k = 5/3. We 
also discuss under which circumstances jets will form internal shock waves for smaller 
deflection angles. 
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1 INTRODUCTION 



In a previous article (Mendoza & Longair 2001), we dis 



cussed a mechanism by which galactic and extragalactic jets, 
can change their original straight trajectory if they pass 
through a stratified cold and high density region. For the 
case of galactic jets this could be a cloud in the vicinity of 
an H-H object. For extragalactic jets, this region could be a 
nearby galaxy, the interstellar medium of the host galaxy or 
the intracluster medium itself. Physically, what is important 
is that the medium interacting with the expanding jet from 
the source has a non-uniform density. 



As mentioned by [eke (1991) and Mendoza & Longair 



(2001 ), when a jet bends, it is in direct contact with the sur- 



roundings and one should expect that entrainment from the 
surrounding gas might cause a severe disruption of the jet 
itself. Assuming that this entrainment is not important, for 
example by an efficient cooling, what its left is a high Mach 
number flow inside a collimated flow that bends. When a su- 
personic flow bends, the character istics emanating at each 
point of the flow tend to int ersect (Landau & Lifshitz 1987 



Courant & Fricdrichs 1976). Since every hydrodynamical 



quantity has a constant value on a given characteristic line, 
this intersection causes the different physical quantities in 
the flow, such as the density or pressure, to be multivalued. 
This situation cannot occur in nature and a shock wave is 
formed. 

The formation of shock waves inside a jet are potentially 
dangerous. These shock waves could give rise to subsonic 
flow in the jet and collimation might no longer be achieved. 
The present article discusses the circumstances under which 
an internal shock wave should be expected in a bent jet. The 



analysis pr esented in this article generalises that made by 
Icke (1991) by introducing relativistic effects into the flow. 
This generalisation has a drastic effect on the results. Rel- 
ativistic jets cannot bend as much as non-relativistic jets. 
As we will see below this occurs because, when relativis- 
tic effects are taken into account, the characteristic lines in 
the flow are beamed in the direction of the flow velocity. 
The beaming increases as the velocity of light is approached 
by the flow. In other words, the chances for an intersection 
between characteristic lines of uniform plane-parallel flow 
(before the bending) and a curved flow increase because of 
this beaming. 

In order to analyse in detail the bending of relativis- 
tic jets we discuss the following points in each section. We 
write down the equations of relativistic hydrodynamics in 
section |^ and their non-relativistic counterparts. This sets 
the scene for discussing the propagation of disturbances in 
flows which leads naturally to the definition of the rela- 
tivistic Mach number and characteristics. This generalises 
the traditional approach to gas dynamics and characteristic 
surfaces su ch as that discussed for n on-relativistic hydrody- 
namics by Landau fc Lifshitz (1987). The r elativistic Mach 
number was first introduced by |Chiu (19751) an d some of its 
properties are well described by IKonigl (1980 ) . The most im- 



portant result to be proved in this section is the beaming of 
characteristic lines in the direction of motion of the flow for 
relativistic flows. Section ^ analyses the relativistic and non- 
relativistic cases for a flow depending on one angular vari- 
able, known as Prandtl-Meyer flow. This leads naturally to 
the generalisation of a type of flow called rarefaction wave s 
in non-relativistic hydrodynamics ( Landau & Lifshitz 1987 ) 
which in turn is essential for the understanding of flows that 
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move through a certain angle. In section |5j we use the main 
results of the Prandtl-Meyer flow to describe steady simple 
waves. These waves form naturally for steady plane parallel 
flow at infinity when the jet turns through an angle with a 
certain curved profile. It is then possible to apply the neces- 



sary physical ingredients to the steady simple waves formed. 



We then calculate an upper limit to the deflection of a jet 
in order to avoid the formation of a shock at the end of its 
curvature. This limit does not mean that a jet which curves 
through a small deflection angle is safe from generating an 
internal shock wave. So, we also calculate a lower limit, for 
which a shock could form at the onset of the curvature. Fi- 
nally, in section [r| we discuss the astrophysical consequences 
of the results described above. 



2 BASIC EQUATIONS 

In order to understand the formation of internal shock waves 
in a bent jet it is necessary to understand some of the basic 
properties of the gas dynamics of relativistic flows. 

The equations of motion for an ideal relativis tic flow arc 
described by the 4-dimensional Euler's equation (Landau & 
Lifshitz 1987|): 



; 8Ui 

dx k 



dp 

dx 1 



: 9p 

dx k 



(1) 



in which Latin indices take the values 0, 1, 2, 3. The vec- 
tor x k — (ct,r), where r represents the three dimensional 
radius vector, c the speed of light and t the time coordi- 
nate. The Galilean metric gik for flat space time is given 
by 500 = 1, gn = gii = 533 = —1 and g ik = when i / k 
The pressure is represented by p, u) = e + p is the enthalpy 
per unit proper volume and e is the internal energy den- 
sity. The four-velocity u k = da; fe /ds where the relativistic 
interval ds is given by ds 2 = gti<ix k dx 1 . The values of the 
different thermodynamic quantities are measured in their 
local proper frame. 

The space components of eq.(^) give the relativistic Eu- 
ler equation: 



dv 

m 



+ v ■ grad v 



-gradp 



v dp 



(2) 



where v is the three dimensional flow velocity and 7 is the 
Lorentz factor for a flow with this velocity. 

For an ideal flow, in the absence of sources and sinks, 



the relativist ic continuity equation is given by (Landau & 
Lifshitz 1987|): 



dx k 



0, 



(3) 



where the particle flux 4-vector n = nu and the scalar n 
is the proper number density of particles in the fluid. 
A polytropic gas obeys the relation: 



p oc n , 



(4) 



where the polytropic index k is a constant and has the value 
5/3 for an adiabatic monoatomic gas in which relativistic 



effects are not taken into account. In the case of an ultrarel- 
ativistic photon gas it has a value of 4/3. It is not difficult to 
show that for a polytropic gas, the speed of sound a and the 
enthalpy per unit mass, the specific enthalp y, w of are re- 
lated to each other by the following formula ( Stanyuokovich 



196C) 



W K - 



1 C 2 



(•») 



The quantities in the relativistic case are defined with 
respect to the proper system of reference of the fluid, 
whereas in classical mechanics these quantities are referred 
to the laboratory frame. In the relativistic case the ther- 
modynamic quantities, such as the internal energy density 
e, the entropy density a and the enthalpy density u are all 
defined with respect to the proper volume of the fluid. In 
non-relativistic fluid dynamics, these quantities are defined 
in units of the mass of the fluid element they refer to. For 
instance, the specific internal energy e, the specific entropy 
s and the specific enthalpy w are all measured per unit mass 
in the laboratory frame. When taking the limit in which the 
speed of light c tends to infinity we must also bear in mind 
that the internal energy density e includes the rest energy 
density nmc 2 , where m is the rest mass of the particular 
fluid element under consideration. Therefore the following 
non-relativistic limits should be taken in passing from rela- 
tivistic to non-relativistic fluid dynamics: 



nn > p\J 1 — v 2 /c 2 S3 p (l — v' 2 /2c 2 ) , 

c — >oo 

2 2 1 2 

e > nmc + pe ~ pc — -pv + pe, 

c — >oo 2 
UJ 1 ( P\ I 2 \ 

> mc + m I e H — ~m(c + w) 

n c->oo \ p J ' 



(6) 



where p is the mass density of the fluid in the laboratory 
frame. 



3 CHARACTERISTICS AND MACH NUMBER 

The properties of subsonic and supersonic flow, are com- 
pletely different in nature. To begin with, let us see how 
perturbations with small amplitudes are propagated along 
the flow for both subsonic and supersonic flows. For simplic- 
ity in the following discussion we consider two dimensional 
flow only. The relations obtained below are easily obtained 
for the general case of three dimensions. 

If a gas in steady motion receives a small perturbation, 
this propagates through the gas with the velocity of sound 
relative to the flow itself. In another system of reference, 
the laboratory frame, in which the velocity of the flow is v 
along the x axis, the perturbation travels with an observed 
velocity u whose x and y components are given by: 



a cos 9 + v 
1 + av cos 6/c 2 

7~ 1 a sin 8 
1 + av cos 6*/c 2 



(7) 
(8) 



according to t he rule for the addition of velocities in spe- 
cial relativity (Landau fc Lifshitz 1994). The polar angle 9 
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Figure 1. Region of influence of small amplitude perturbations. 
A perturbation of small amplitude is produced in the flow at some 
point 0. This is carried by the flow which has velocity v. In the 
case of subsonic flow, as shown in panel (a), the perturbation is 
able to propagate to the whole space. When the flow is supersonic 
the perturbation is propagated only downstream inside a cone 
with aperture angle 2a. The speed of sound a and the angle 8 
are measured in a frame of reference in which the flow is at rest 
-the proper frame of the flow. The vector u is the relativistic 
addition of vectors v and ae! r , where e! r is a unit radial vector in 
the proper frame of the flow. 



and the velocity of sound a are both measured in the proper 
frame of the fluid. Since a small disturbance in the flow 
moves with the velocity of sound in all directions, the pa- 
rameter 9 can have values 9 ^ 2n. This is illustrated 
pictorially in Fig. [lj. 

Let us consider first the case in which the flow is sub- 
sonic, as illustrated in case (a) of Fig. Since_by_defini- 



it follows from eqs.(^)-(^) that 
- 0. In other words, the region 



tion v < a and c > av. 
u x {9 = 7r) < and u y (9 = i 
influenced by the perturbation contains the velocity vector 
v. This means that the perturbation originating at is able 
to be transmitted to all parts of the flow. 

When the velocity of the flow is supersonic, the situa- 
tion is quite different, as shown in case (b) of Fig. [I]. For 
this case, u x (9 = tt) > and u y (9 =7r) = 0. In other words, 
the velocity vector v is not fully contained inside the region 
of influence produced by the perturbation. This implies that 
only a bounded region of space will be influenced by the per- 
turbation originating at position 0. For the case of steady 
flow, this region is evidently a cone. Thus, a disturbance 
arising at any point in supersonic flow is propagated only 
downstream inside a cone of aperture angle 2a. By defini- 
tion, the angle a is such that it is the angle subtended by 
the unit radius vector e r with the velocity vector v at the 
point in which the azimuthal unit vector e a is orthogonal 
with the tangent vector d(ae' r )/d9 to the boundary of the 
region influenced by the perturbation. The unit vector e' r 
is the unit radial vector in the proper frame of the flow. In 
other words, the angle a obeys the following mathematical 
relation: 



On the other hand, since tana 
eqs.(|)-(|) and eq.@ that: 



u y /u x , it follows from 



tan( 



1 - 



and so eq.fllOD gives a relation between the angle a, the 
velocity of the flow v and its sound speed a: 



tan a = 7 



a/v 



(11) 



This variation of the angle a with the velocity of the 
flow is plotted in Fig. ^ for the case in which the gas is 
assumed to have a relativistic equation of state, that is, when 
p — e/3. The important feature to note from the plot is that 
the aperture angle of the cone of influence is reduced when 
the velocity of the flow approaches that of light. 

From eq.(pr|) it follows that, as the velocity of the flow 
approaches that of light, the angle a vanishes. In other 
words, as the velocity reaches its maximum possible value, 
the perturbation is communicated to a very narrow region 
along the velocity of the flow. 

In studies of supersonic motion in fluid mechanics it 
is very useful to introduce the dimensionless quantity M 
defined as: 



1 "la a 
— = sin a = , 

M j v 



(12) 



according to eq.(|Tl]). The quantity j a = l/y/T— (a/c) 2 is the 
Lorentz factor calculated with the velocity of sound a. The 
number M has the property that M — > 1 as v — > a and M — > 
oo as v — > c. It also follows that M > 1 if and only if v>a. 

The surface bounding the region reached by a distur- 
bance st arting from the origin is called a characteristic 
surface ( Landau fc Lifshitz 1987 ) . In the general case of ar- 
bitrary steady flow, the characteristic surface is no longer 
a cone. However, exactly as it was shown above, the char- 
acteristic surface cuts the streamlines at any point at the 
angle a. 

Let us briefly discuss the non-relativistic limit of the 
different physical circumstances presented above. To do this, 
we use the relations in eq.(^|) with c^oo and, as usual for 
the non-relativistic case, we represent the speed of sound by 
c. 

The dimensionless number M satisfies the following re- 
lation: 



M 



(13) 



and is called in non-relativistic hydrodynamics the Mach 
number^ 



d9 Sa ~ U ' 
Substitution of eqs.(Q)-(^|) into eq.(^|) gives: 



tana 



1 av 



tan t 



0) 



(10) 



1 The relativistic generalisation of the Mach n umber as presented 
in eq. (|l2|) was first calculated by |Chiu (1973 ), who reduced the 



problem of steady relativistic gas dynamics to an equivalent non- 
relativistic flow. From eq. (jHj) it follows that this number is in fact 
a definition of the proper Mach number since it is defined as the 
ratio of the space component of the relativistic four-velocity of 
the flow jv to the same com nonent of the relativistic four-velocity 
of sound 7 a a (Konigl 1980). 
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Figure 2. Region of influence of a perturbation for different 
values of the velocity of a relativistic gas with a sound speed 
a = c/\/3. From left to right the closed loops correspond to values 
of the velocity v of 0.0, 0.2, 0.4, . . . , 0.8 in units of the speed of 
light c. The perturbation is assumed to originate at the origin of 
the proper system of reference of the flow. In the case of super- 
sonic flow, the region of influence occurs only downstream inside 
a cone with aperture angle 2a. This cone surrounds the corre- 
sponding loop and is tangent to it. When the flow is subsonic, 
the perturbation is transmitted to all the flow. As a particular 
example, the cone and the velocity vector have been drawn for 
the case in which the velocity is u = 0.8c. 



The results obtained concerning the relativistic and 
non-relativistic Mach number M can be rewritten in the 
following way: the dimensionless Mach number M increases 
without limits as the velocity of the flow takes its maximum 
possible value. This maximum value is the speed of light 
in the relativistic case and infinity in the non-relativistic 
one. The Mach number tends to zero as the velocity of the 
flow vanishes, and tends to unity as the velocity of the flow 
tends to the velocity of sound. The Mach number is greater 
than one for supersonic flow and less than unity when the 
velocity of the flow is subsonic in both, the relativistic and 
non-relativistic cases. 



dv r 



v<t> 



, d(u07i7) 
uiy/n — const, 



(14) 

(15) 
(16) 



where v r and v$ are the components of the velocity in the 
radial and azimuthal direction s respectively. Eq.(|l6[) is the 
relativistic Bernoulli equation (Landau & Lifshitz 1987) for 
this problem. 

Using the definition of the speed of sound, 



dp 

7k: 

eq.(|l5|) and eq.(|l6|) it is found that: 



v r + 



+ iv 1 



a 2 d(j> 



ln(w/n) = 0. 



(17) 



(18) 



On the other hand, differentiation of 7 2 with respect to the 
azimuthal angle </> and using eq.(^), gives: 



, , „ + *?±) + e [ ! _ v l±<) d ln (w/n) = a (19) 



Multiplication of eq.(|lq) by and substracting this from 



ultrp 
eq.(|ig|) gives 



■2 1 1 . Vt 
= a I 1 



(20) 



Bernoulli's equation, eq.(|16[), together with the value of 
the specific enthalpy for a polytropic gas given in eq. (jHj) , can 
be rewritten 



K-1C 2 



= 1- 



1 »r 2 



K-1C 2 



(21) 

in which it has been assumed that at some definite point, 
the flow velocity vanishes and the speed of sound has a value 
ao there. It is always possible to make the velocity zero at a 
certain point by a suitable choice of the system of reference. 
Eqs.(po|)-(|2l|) can be solved in terms of v r and v^: 



4 PRANDTL-MEYER FLOW 

Let us describe briefly the exact solution of the equations 
of hydrodynamics for plane steady flow which depends only 
on the angular variable <f> only. This proble m was first in- 



vestigated by Prandtl and Meyer in 1908 ( Landau fc Lif- 
shitz 1987; Courant fc Friedrichs 1976) for the case in which 



relativistic effects were not taken into account. The full rel- 



ativistic solution to the problem is due to Kolosnitsyn & 
Stanyukovich (1984j ). 

For this case, Euler's equation, eq. (Jl|) and the continuity 
equation, eq.tel) can be written: 



where 



F\a) = 1 



v'r/c 2 = 1 — F 2 (a), 
v\ = a 2 F 2 (a), 



-I 2\ 2 / 2\ "I 



K-1C 2 



1 a' 



(22) 
(23) 



2\ "2 



(24) 

Because v r dy r = c 2 F(a)F' (a)da, cq.(|l^) give s the re- 
quired solution (Kolosnitsyn & Stanyukovich 1984): 



■ + 4>o — ±c 



F'{a)da 



(25) 
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This equation gives the speed of sound as a function of the 
azimuthal angle. From eqs.(p2[)-(p^) it follows that the radial 
and azimuthal velocities can be obtained as a function of the 
same angle <f>. As a result, all the remaining hydrodynamical 
variables can be found. The sign in eq.fcq) can be chosen 
to be negative by measuring the angle <f> in the appropriate 
direction and we will do that in what follows. 

Let us consider now the case of an ultrarelativistic gas 
and integrate eq.(M) by parts, to obtain: 




' + 4>o = — arccos F(a) + c / ^ arccos F(a). (26) 
a J a 1 

For the case of an ultrarelativistic gas, the speed of sound 

\fk~- 



a is given by ( Btanyuokovich 1960 ) a = <Jk — 1 c . In other 
words, this velocity is constant and so the integral in eq.(|26[) 
is a Lebesgue integral. Since this integral is taken over a 
bounded and measurable function over a set of measure zero, 
its value is zero. 

Using eqs.(p^~)-(p3[) and eq.(p6f) the des i red solution is 
obtained ( Kolosnitsyn fc Stanyukovich 1984 ; Konigl 1980|) : 



v r — c sm 



{Vk- 1(0 + o )} 



V</> 



Vft — 1 c cos {vk - 1 (4> + 0o)} 



(27) 
(28) 



for an ultrarelativistic equation of state of the gas. 

For the non-relativistic case, in which c — > oo, eq.(|25| 
gives for a polytropic gas with polytropic index k: 



■ + 00 



K+l 



K-1J ^ C 5 



where 



CV K-l 



K-1C 2 



and so, t he required solution is ( Kolosnitsyn fc Stanyukovich 
1984)71 



K + l fc 
■ + 4>o = \ 7 arccos — 

K — 1 V C* 



(29) 



where the gpppH nf gnnnrl n Viaa hppti rpwrittpn ag r tn hp 



consistent in the non-relativistic case. The critical velocity 



of sound c* is given by ( Landau & Lifshitz 1987 



2 



2 

-Co- 



(30) 

The value for the velocities can thus be calculated from 
eq.(0) and eq.@ with F(a) = l: 



K+l 



K-l 
K+l 



c* cos 



K- 1 



(31) 



(32) 



K+l 

Some important inequalities must be satisfied for the 
flow under consideration. First of all, eq.(^i|) together with 
eq.(^) and the first law of thermodynamics , d(uj/n) = 
TA(a/n) + (l/n)dp ( [Landau fc Lifshitz 1987| ), imply that 
dp/d(f> < 0. Using this inequality and the fact that de = 



Figure 3. Relation between the velocity vector v = v r e r +ti^e0 
and the angle \> as a function of the azimuthal angle <f>. \ is the 
angle that the velocity vector makes with certain fixed axis with 
origin O. 



c 2 dp/a 2 combined with the first law of thermodynamics, it 
follows that dn/dif) < 0. Also, using eqs.(|23)-(E3) it follows 
that dv/d(f> oc — da/d<f> and necessarily dvjd(f)>0. 

On the other hand, the angle % tna t the velocity vector 
makes with some fixed axis is related to the velocity and the 
azimuthal angle <f> by: 

X = (j> + arctan (w^ /v r ) (33) 

as it shown in Fig. ^[ Thus, since the <f> component of Euler's 
equation, eq.(|l|) implies that: 



V r + 



d(f> J n dn 



= 0, 



it follows that dx/d(j> — — (y 2 ^ 2 ui/c 2 ) _1 dp/d(j>. 

In other words, we have proved that for the flow with 
which we are concerned, the following inequalities are satis- 
fied: 

dp/d(l> < 0, dn/d<j) < 0, dv/dcf) > 0, dx/d(j) > 0. (34) 

A flow with these properties is describ ed as a rarefaction 
wave in non-relativistic fluid dynamics ( Landau & Lifshitz 



1987) and we will use this name in what follows. 

Another, very important property of this rarefaction 
wave is that the lines at constant (j> intersect the stream- 
lines at the Mach angle, that is, they are characteristics. 
Indeed, from Fig. ^, it follows that the angle a between 
the line (j> = const and the velocity vector v is given by 
sin a = Vcfy/v. Using eqs.(p2"l)-(p4j) it follows that this rela- 
tion can be written as eq.(Q2). Because all quantities in the 
problem are functions of a single variable, the angle cj>, it fol- 
lows that every hydrodynamical quantity is constant along 
the characteristics. 



5 STEADY SIMPLE WAVES 

Let us now consider the two dimensional problem of steady 
plane parallel flow which then turns through an angle as 
it flows round a curved profile. A particular case of this 
problem, when the flow turns through an angle is described 
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by Landau & Lifshitz (1987). For this particular situation 



the Prandtl-Meyer flow is obviously the solution and so, the 
hydrodynamical quantities depend on a single variable, the 
angle (j> measured from a defined axis at the onset of the 
curvature. Because of this, all quantities can be expressed 
as functions of each other. Since this case is a particular 
solution to the general problem, it is natural to seek the 
solutions of the equations of motion in which the quanti- 
ties p, n, v x , v y can be expressed as a function of each 
other. Evidently this imposes a restriction on the solution of 
the equations of motion since for two dimensional flow, any 
quantity depends on two coordinates, x and y, and so any 
chosen hydrodynamical variable can be written as a function 
of any other two. 

Because of the fact that the flow is uniform at infinity, 
where all quantities are constant, particularly the entropy, 
and because the flow is steady, the entropy is constant along 
a streamline. Thus, if there are no shock waves in the flow, 
the entropy remains constant along the whole trajectory of 
the flow and in what follows we will use this result. 

In this case, Euler's equation, eq.(^), and the continuity 
equation, eq.feh, are respectively: 




Figure 4. Supersonic uniform flow at the left of the diagram 
bends around a curved profile OH. The Mach angle ol\ is the 
angle made by the characteristics and the streamlines before the 
onset of the curvature. The characteristics make an angle </>i with 
the "continuation" of the rarefaction wave formed at the onset of 
the curvature and the angle cf> is measured from the line OA'. The 
curvature causes the characteristic lines to intersect eventually 
and this occurs at point K in the diagram, giving rise to a shock 
wave represented as the segment KL. 



dv x , dv, 

Vx — h Vy 



dx 



dy 



8Vy 8Vy 

Vx-tT+Vy-^- 

dx ay 

JL (w) + JL (wi) 



c 2 dp 

7W dx ' 

e_dp_ 

~fij dy ' 



Rewriting these equations as Jacobians^] we obtain: 



d(vx,y) „, d(v x ,x) 
' d(x,y) 

d(y y ,y) 

d(x,y) 



d(x,y) 
d(v y ,x) 
V d(x,y) 



= + 



c 2 d(p,y) 
juj d(x,y)' 

c 2 d{p,x) 
-yuj d(x,y) ' 



d{*iv x n,y) _ d(^v y n,x) 
d(x,y) d(x,y) 



: 0. 



We now take the coordinate x and the pressure p as in- 
dependent variables. To make this transformation we have 
to multiply the previous set of equations by d(x,y)/d(x,p). 
This multiplication leaves the equations the same, but with 
the substitution d(x,y) — > d(x,p). Expanding this last re- 
lation and because all quantities are now functions of the 
pressure p but not of a;, it follows that: 



dy\ dv x 



dx J dp 
dy\ dvy ; 
dx J dp 



dx ) dp 



c?_dy 

7u> dx ' 

c 2 dy 

yui dx ' 



dVy_ 

dp 



dy dv : 
dx dp 



Here we have taken dy/dx to mean the derivative 



The Jacobian d(a,b)/d(x, y) is defined as: 



8(a, b) 
d(x,y) 



det 



da/dx 
db/dx 



da/dy 
db/dy 



at constant pressure: (dy/dx) p . Since every hydrodynamic 
quantity is assumed to be a function of the pressure, then 
in the previous set of equations it necessarily follows that 
dy/dx is a function which depends only on the pressure, 
that is (dy/dx) p = fi{p)- Therefore: 



V = xfi(p) + hip)- 



(35) 



No further calculations are needed if we use the solution 
for the case in which a rarefaction wave is formed when flow 



turns through an angle ( Landau fc Lifshitz 1987 ). This so- 
lution is given by the results of section ^. As was mentioned 
in that section, all hydrodynamical quantities are constant 
along the characteristic lines (f> — const. The particular so- 
lution of the flow past an angle obviously corresponds to 
the case in which / 2 (p)=0 in eq.(^). The function A (p) is 
determined from the equations obtained in section H. 

For a given constant value of the pressure p, eq.(^), 
gives a set of straight lines in the x-y plane. These lines 
intersect the streamlines at the Mach angle. This occurs be- 
cause the lines y — xfi (p) for the particular solution of the 
flow through an angle have this property. In other words, 
one family of characteristic surfaces correspond to a set of 
straight lines along which all quantities remain constant. 
However, for the general case, these lines are no longer con- 
current. 

The properties of the flow described above are anal- 
ogous t o the non-relativistic eq uivalent known as simple 
waves (Landau & Lifshitz 1987). In what follows we will 
use this name to refer to such a flow. 

Let us now construct the solution for a simple wave once 
the curved profile is fixed. Consider the profile as shown in 
Fig. ^. Plane parallel steady flow streams in from the left of 
point O and flows around the curved profile. Since we as- 
sume that the flow is supersonic, the effect of the curvature 
starting at O is communicated to the flow only downstream 
of the characteristic OA generated at point O. The char- 
acteristics to the left of OA, region 1, are all parallel and 
intersect the x axis at the Mach angle ai given by eq.(^2|): 
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sm«i = 



yjl- (vi/c) 2 a 
\A - (o/c) 



2 Wl 



(36) 



where the velocity vi is the velocity of the flow to the left 
of the characteristic OA. In eqs,([25[)-(p^) the angle of the 
characteristics is measured with respect to some straight 
line in the x-y plane. As a result, we can choose for those 
equations the constant of integration 0o = 0. This means 
that the line from which the angle is measured has been 
chosen in a rather special way. In order to find the line which 
is the characteristic for — 0, let us proceed as follows. 
When = and the gas is ultrarelativistic, eqs.(^7|)-(|28|) 
show t hat the velocity v = a and for the non-rclativistic 



case, it follows from eqs.'f^ )-([j2[) that the velocity takes the 
value v = c. In both cases this means that the line — 
corresponds to the point at which the flow has reached the 
value of the local velocity of sound. This, however, is not 
possible since we are assuming that the flow is supersonic 
everywhere. Nevertheless, if the rarefaction wave is assumed 
to extend formally into the region to the left of OA, we 
can use these relations and then the characteristic line must 
correspond to a value of given by: 



K + 1 Ci 

arccos — 

K — 1 V c* 



for a non-relativistic gas according to eq.(E9|), and 



(37) 



(38) 



c Vl-(Vc) 2 
si = — arccos — — . 

v/i - («/c) 2 ' 



for the ultrarelativistic case according to eq.(|21[), eq.(|24[) 
and eq.(^). The angle between the characteristics and 
the x axis is then given by: 0* — 0, where 0* = cti + 
0i, and the angle a± is the Mach angle in region 1. The 
x and y velocity components in terms of the azimuthal angle 
9 are given by: 



v x = v cos f 



v v = v sm ( 



(39) 



and the values for the magnitude of the velocity, the angle 
9 and the pressure are given by: 



1 + - 



K + l 



arctan 



ft + 1 



■ cot 



K+l 



p* COS 



2k/(k-1) 



K+l 



0, 



(40) 



(41) 



(42) 



for a non-relativistic gas according to eqs.(|29|)-(p2|) and us- 
ing the fact that the Poisson adiabatic for a polytropic gas 
means that: pc~ 2K '^ n ~^ =const. In the case of an ultrarela- 
tivistic gas, eqs.(|2^)-(p8[) together with Bernoulli's equation 
and the fact that the enthalpy density w = /tp/ (k — 1) give: 



v 2 — c 2 { 1 — (2 — k) cos 2 \fk — 10} 



arctan • 



I k — \ cot V 'ft — 10} , 

/o \ — re/2(re— 11 —k/(k— 11 / T, 

p = po (2 — k) ' * 'cos /v 'y/K—ltp. 



(43) 
(44) 
(45) 



Since the angle 0* — is the angle between the char- 
acteristics and the x axis, it follows that the line describing 
the characteristics is: 



y — a; tan (0* — 0) + G(0). 



(46) 



The function G(0) is obtained from the following argu 



merits for a given profile of the curvature ( Landau & Lifshitz 



1987). If the equation describing the shape of the profile is 
given by the points X and Y where Y = Y(X), the velocity 
of the gas is tangential to this surface, and so: 



tanf? ; 



AY 
AX' 



(47) 



Now, the equation of the line through the point (X, Y) which 
makes an angle 0* — with the x axis is: 

y - Y = (x - X) tan (0» - 0) . (48) 
Eq.(|48[) is the same as eq.(^) if we set: 

G(0) = Y- Atan(0, -0). (49) 

If we start from a given profile Y = Y(X) then, using 
eq.(^) we can find the parametric set of equations: X — 
X(9), and Y = Y(9). Substitution of = 0(0) from eq.@ or 
eq.(|44]) depending on whether the gas is non-relativistic or 
ultrarelativistic, we find X = X((j>) and Y = Y(<f>). Substitu- 
tion of this in eq.(^) gives the required function G(0). 

If the shape of the surface around which the gas flows 
is convex, the angle 9 that the velocity vector makes with 
the x axis decreases downstream. The angle — 0* between 
the characteristics leaving the surface and the x axis also 
decreases monotonically. In other words, characteristics for 
this kind of flow do not intersect resulting in a continuous 
and rarefied flow. 

On the other hand, if the shape of the surface is con- 
cave as shown in Fig. m, the angle 9 increases monotoni- 
cally and so does the angle the characteristics make with 
the x axis. This means that there must exist a region in 
the flow in which characteristics intersect. The value of the 
hydrodynamical quantities is constant for every character- 
istic line. This constant however changes for different non- 
parallel characteristics. In other words, at the point of in- 
tersection different hydrodynamical quantities, for example, 
the pressure, are multivalued. This situation cannot occur 
and results in the formation of a shock wave. This shock 
wave cannot be calculated from the above considerations, 
since they were based on the assumption that the flow had no 
discontinuities at all because the entropy was assumed to be 
constant. However, the point at which the shock wave starts, 
that is point K in Fig. W, can be calculated from the follow- 
ing considerations. We can work out the inclination of the 
characteristics as a function of the coordinates x and y. 
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This function <j>(x, y) becomes multivalued when these coor- 
dinates exceed certain fixed values, say xo and yo- At a fixed 
x = xq the curve giving the value of as a function of y be- 
comes multivalued. That is, the derivative (d(j>/dy) x =00, 
or (dy/d(j)) x = 0. It is evident that at the point y = yo the 
curve 4>{y) must lie in both sides of the vertical tangent, else 
the function <j>(y) would already be multivalued. This means 
that the point (xo, yo) cannot be a maximum, or a minimum 
of the function <f>(y) but it has to be an inflection point. In 
other words, the coordinates of po int K in Fig. H can be cal- 



culated from the set of equations ( Landau fc Lifshitz 1987 ) 



where 



^/(k — l)/(/t + 1) if the gas is non-relativistic, 



for an ultrarelativistic gas. 



(52) 



As was mentioned above, if the jet is sufficiently nar- 
row, it appears that it can safely avoid the formation of an 
internal shock. However, differentiation of eq.(^) with re- 
spect to the angle the velocity vector makes with the x axis, 
that is the deflection angle 9, implies that: 



dy 
do 



= 0, 



\d<l> 2 



= 0. 



(50) 



When the profile is concave, the streamlines that pass 
above the point O in Fig. ^ pass through a shock wave and 
the simple wave no longer exist. Streamlines that pass be- 
low this point seem to be safe from destruction. However, 
the perturbing effect from the shock wave KL influences this 
region also, and so it is not possible to describe the flow there 
as a simple wave. Nevertheless, since the flow is supersonic, 
the perturbing effect of the shock wave is only communi- 
cated downstream. This means that the region to the left 
of the characteristic PK, which corresponds to the other set 
of characteristics emanating from point P, does not notice 
the presence of the shock wave. In other words, the solution 
mentioned above, in which a simple wave is formed around a 
concave profile is only valid to the left of the segment PKL. 



6 CURVED JETS 

Let us now use the results obtained in sections ^ and ^ and 
apply them to the case of jets that are curved due to any 
mechanism, for example the interaction of the jet with a 
cloud as was d iscussed by us in a previous paper ( Mendoza & 
Longair 2001 ), or due to the ram pressure of the intergalactic 



gas as in the core of radio trail sources. 

The greatest danger occurs when the jet forms inter- 
nal shock waves. This is because, after a shock, the normal 
velocity component of the flow to the surface of the shock 
becomes subsonic and the jet flares outward. Nevertheless, 
as we have seen in section ^, the shock that forms when gas 
flows around a curved profile (such as a bent jet due to ex- 
ternal pressure gradients) does not start from the boundary 
of the jet. It actually forms at an intermediate point to the 
flow. In other words, it is possible that, if a jet does not bend 
too much the intersection of the characteristic lines actually 
occurs outside the jet and the flow can curve without the 
production of internal shocks. 

As we have seen in section ^ the Mach angle of the flow, 
relativistic and non-relativistic, does not remain constant in 
the bend (see for example eq.(^) and eq.(^)). The Mach 
number monotonically decreases as the bend proceeds. 

Eq.dO) and eq.@ imply that: 



do 
AO 



1 + 



r + 1 

M 2 - 1 



with 



r = 



k if the gas is non-relativistic, 

k/(2 — k) for an ultrarelativistic gas. 



(53) 



(54) 



The Mach number M is given by eq.(]l3|) and eq.(|l^) re- 
spectively. As the Mach number M — > 1, then the deriva- 
tive Aa/A9~ >oo. This means that the rate of change of the 
Mach angle with respect to the deflection angle grows with- 
out limit as the Mach number decreases and reaches unity. 
On a bend, the Mach number decreases and care is needed, 
or else the characteristics will intersect at the end of the cur- 
vature. There is only one special shape for which this effect 
is bypassed and this occurs wh en the increase of 9 matche s 
exactly with the increase of a (Courant & Friedrichs 1976), 



but of course, this is quite a unique case. It appears however, 
that whatever the thickness of the jet it cannot be bent more 
than the point at which Aa/A9 exceeds the rate of change 
of 9 with respect to the bendin g angle 9. In other words, 
Aa/A9 <A9/A9 = 1 ( tcke 199l|) . From this last inequality 
and eq.fca) a value of the Mach number can be obtained: 



M* 



V3- 



(55) 



If the Mach number in the jet decreases in such a way that 
the value M* is reached, then a terminal shock is produced 
and the jet structure is likely to be disrupted. It is important 
to note that this terminal shock is weak since M> 1 and so, 
it might not be too disruptive. Nevertheless, this monotonic 
decrease of the Mach number makes the jet flare outwards, 
even if the terminal shock is weak. 

Let us now calculate an upper limit for the maximum 
deflection angle for which jets do not produce terminal 
shocks. In order to do so, we rewrite eq.(pl|) in the following 
way: 



• 1 1 

—9 = arcsm -— H arctan 

M n 



{ M \/A/ 2 -l}-^ 



(56) 



tan a = —ft cot fj, (a + 1 



(51) 



To eliminate the constant 0* from all our relations, we can 
compare the angle 9 evaluated at the minimum possible 
value of the Mach angle M = M+ with 9 evaluated at its 
maximum value M = 00. In other words, the angle # ma x 
defined as: 
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Figure 5. Sketch of a curved jet of radius D that develops a 
shock at the beginning of the curvature. The curve is assumed to 
be a circle with radius R, so that this approximation is valid on 
a sufficiently small region about the onset of the curvature. The 
Mach angle of the jet is a at the left of the characteristic OA that 
emanates from the point where the bending starts. 



(M = M*)-0(M=oo) 

74.21 ° for the non-relativistic case, 
47.94 ° for an ultrarelativistic gas. 



(57) 



are upper limits to the deflection angle. Jets which bend 
more than this limiting value max develop a terminal shock. 

This upper limit however, does not mean that the jet 
is immune from developing an internal shock if it is bent by 
a smaller angle. Indeed, let us suppose that the jet bends 
and that the curvature it follows is a segment of a circle as 
it is shown in Fig. ^. The circle can be considered to be the 
circle of curvature of the jet's trajectory formed at the onset 
of the bending. According to the figure, the equation of the 
characteristic OA that emanates from the point O, where 
the curvature starts is: 



y — x tan a 



(58) 



Once the flow has curved d8 degrees, the characteristic 
at this point is given by: 

y = (x — Rdd) tan(a + da + d<f>) 
~ x tan a + x (da + dO) / cos 2 a — Rd6 tan a, 

where R is the radius of curvature of the circular trajec- 
tory. The intersection of this characteristic and that given 
by eq.(S) occurs when the y coordinate has a value: 



D = 



i?sin a 
1 + da/d9' 




Substitution of eq.feSh gives ([eke 1991) 



15 20 25 30 35 40 45 50 



Figure 6. Plot of the maximum ratio D/R as a function of the 
difference 9 — 9+ where 9 is the deflection angle and 9+ is the 
maximum bending angle a jet can have in order not to produce a 
terminal shock. The plot refers to the points for which a shock at 
the beginning of the curvature (which was assumed to be a circle) 
has developed. Jets with parameters which lie below the curve in 
any case do not develop any internal shocks at all for this partic- 
ular circular trajectory. The plot at the top was calculated using 
the results in which the gas is non-relativistic and its polytropic 
index is 5/3. The plot at the bottom was made by considering 
the gas to be ultrarelativistic and relativistic effects in the bulk 
motion of the flow were taken into account. For this second plot, 
the polytropic index was assumed to be 4/3. The numbers in each 
plot correspond to the values of the Mach number in the flow. 



between the deflection angle 8 and the maximum deflection 
angle 8+ = 8(M+), as is shown in Fig. [| 

Jets for which the ratio D/R lies below the curve do not 
develop any shocks at all. For example, consider a jet with 
a given Mach number for which its ratio D/R is given. As 
the width of the jet increases (or the radius of curvature of 
the profile decreases), it comes a point in which a shock at 
the onset of the curvature is produced. In the same way, jets 
with a fixed ratio D/R for a given Mach number which are 
initially stable, so that they lie below the curve, can develop 
a shock at the beginning of the curvature by increasing the 
bending angle of the curve. 



^ = 7t~ — (A/ 2 - 1) M~ 4 . (59) 

r r+ i 1 ' v ' 

Using eq.(|56|) and eq.(^) it is possible to make a plot 
in which two zones separate the cases for jets which develop 
shocks at the onset of the curvature, and the ones that do 
not. Indeed, we can plot the ratio of the width of the jet D 
to radius of curvature J? as a function of the difference 8 — 8+ 



7 DISCUSSION 

The relativistic Mach angle is smaller for a given value of 
the velocity of the flow than its non-relativistic counterpart 
as was proved in section ^[ This fact is extremely important 
when analysing the possibility of the intersection of different 
characteristics in a bent jet. This intersection is what gives 
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rise to the creation of shock waves. For a relativistic flow, the 
characteristics, which make an angle equal to the Mach angle 
to the streamlines, are always beamed in the direction of 
the flow. Thus, when a jet starts to bend, the possibilities of 
intersection between some characteristic line in the curved 
jet and the ones before it curves, become more probable than 
for their non-relativistic counterpart. 

This difference results in a severe overestimation of the 
maximum bending angle 9 mal when a non- relativistic treat- 



ment is made to the problem. For example, Icke (1991) used 



the non-relativistic analysis in the discussion of the gener- 
ation of internal shocks due to bending of jets. Using the 
non-relativistic equations described above, but with a poly- 
tropic index k = 4/3, the value for the maximum deflection 
angle i s ffima.y — 134.16 . This is much greater than the value 



of # max =74.21 ° obtained with a full relativistic treatment. 



The analysis made by Icke (1991 ) is important for jets in 
which the microscopic motion of the flow inside the jet is rel- 
ativistic, but the bulk motion of the flow is non-relativistic. 



Radio trail sources (Begelman et al. 1984) show consid- 
erable bending of their jets with deflection angles of about 
90 ° in many cases. Since the bending is produced by the 
proper motion of the host galaxy with respect to the in- 
tergalactic medium, the deflection angle cannot be greater 
than 90°. The results presented in eq.(^) show that jets 
which have a relativistic equation of state and a bulk rela- 
tivistic motion of the gas within its jet, cannot be deflected 
more than ~ 50 ° . Since the deflections of radio trail sources 
are greater than this value, this result would imply that 
most radio trail sources should generate shocks at the end 
of their curvature. However, observations (s ee for example 

and references 



Eilek et al. 1984 



O'Dea 1985 



ie Young 1991 



within) show that the velocity of the material of the jets 
< 0.2-0. 3c. Therefore, the bulk motion of the flow is non- 
relativistic, despite the fact that the gas inside the jet has a 
relativistic equation of state. As we saw before, this implies 
that the value of the maximum angle is # max = 134.16 °. 

In other words, these type of jets develop a terminal 
shock if their jets bend more than ~ 135°. This seems to 
be the reason why radio trail sources are able to bend so 
much without resulting in an internal shock wave that could 
potentially cause disrupt ion of its structure. 

In a previous paper (Mendoza & Longair 2001), we dis- 
cussed the possibility of a bent jet in the rad io galaxy 3C 34 
using the observations of Best et al. (1997). According to 
these observations, the radio source lies more or less in the 
plane of the sky, and so, if the western radio jet is curved, 
this must be of the order of 10 ° . The value 6 ~ 10 ° is well 
below the upper limit of ~ 50 ° calculated in eq.(^), so that 
no terminal shock would be produced by the deflection of 
the jet. From the lower plot of Fig. ^| and because the angle 
6* <C 1 for a high relativistic flow, it follows that if the tra- 
jectory of the jet in 3C 34 is circular, then in order not to 
produce an internal shock at the onset of the curvature, the 
ratio D/R has to be less than ~0.08. 

In the analysis made above, we have calculated how 
shocks can be generated inside a bent jet. These shocks are 
special in the sense that they do not reach the surface bound- 
ary of the jet. Instead they are generated away from the 
walls of the jet. For jets with an ultrarelativistic equation 
of state that possess a relativistic bulk motion, a shock is 
internally generated if they bend more than ~ 50 ° . If their 



bulk motion is non-relativistic, the shock is generated when 
the bending angle is more than ~ 135 ° . Jets with a poly- 
tropic index of 5/3 that move non-relativistically generate 
a shock if the bending angle exceeds the value ~ 50 ° . These 
angles are only upper limits and the precise conditions under 
which a shock is produced have to be calculated individu- 
ally. However using the diagrams presented in Fig. and 
that presented by (Icke 1991) it is possible to see if a shock 



is produced at the onset of the curvature. 

All radio sources in which bendings of radio jets have 
been observed appear to satisfy the upper limits discussed 
above. So, it seems that jets are perhaps unstable if an in- 
ternal shock is generated in a curvature. However, various 
observations and theoretical work in galactic and extragalac- 



tic sources (see for example Canto et al. 198S ; Falle fc Raga 



1992 , 1995; Komissarov & Falle 1998, and references within) 
show that internal shocks within a jet are a good mechanism 
that, under certain circumstances, can collimate the jets. 

The question of whether an internal shock will lead to 
disruption of the jet is unknown and still a matter of debate. 
We aim to give to the problem an answer in a future paper. 
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